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Abstract 

We study the issue of duality violations in the VV — AA vacuum polarization 
function in the chiral limit. This is done with the help of a model with an 
expansion in inverse powers of the number of colors, N c , allowing us to consider 
resonances with a finite width. Due to these duality violations, the Operator 
Product Expansion (OPE) and the moments of the spectral function (e.g. the 
Weinberg sum rules) do not match at finite momentum, and we analyze this 
difference in detail. We also perform a comparative study of many of the different 
methods proposed in the literature for the extraction of the OPE parameters and 
find that, when applied to our model, they all fare quite similarly. In fact, the 
model strongly suggests that a significant improvement in precision can only be 
expected after duality violations are included. To this end, we propose a method 
to parameterize these duality violations. The method works quite well for the 
model, and we hope that it may also be useful in future determinations of OPE 
parameters in QCD. 



1 Introduction 

The Operator Product Expansion (OPE) lj is of paramount importance in QCD 
because it is the short-distance window not only to fundamental parameters of the 
Lagrangian such as quark masses and the coupling constant a s , but also to other im- 
portant quantities not explicitly appearing in the Lagrangian, such as the condensates. 
Moreover, there are cases, like that of the Hyv ~ ^-aa vacuum polarization, where 
these condensates turn out to be related to electroweak coupling constants governing 
the physics of kaon decay and its associated CP violation [2] . It is therefore essential to 
understand how the OPE works and what its properties of convergence are. It is clear 
that using the expansion in a region without good convergence (or with no convergence 
at all) may result in errors in the determination of these important parameters. 

The OPE in QCD is believed to be only an asymptotic expansion in inverse powers 
of q 2 in the complex q 2 plane excluding the Minkowski region, Re q 2 > 0. However, it 
is precisely in the Minkowski region where experimental data are available, so that in 
order for the OPE to make contact with the experimental data, one has to learn how 
to make the necessary rotation in the complex plane. A famous trick J3], relying on 
Cauchy's theorem, relates the integral of the spectral function (i.e. of the experimental 
data) pit) = (1/tt) Im Il(t) to that of the corresponding Green's function IT(g 2 ) of 
complex argument q 2 , over the contour of Fig. ^ chosen counter-clockwise. One finds 



where P(x) is any polynomial. As it stands, Eq. is of course an exact mathematical 
statement. The right-hand side of Eq. (0), however, contains the full Green's function 
Il(g 2 ) which is not available in QCD. It is then assumed that, if the contour radius sq 
is large enough, use of the expansion UopE{q 2 ) instead of the exact function will be a 
good enough approximation. This is the rationale underlying all modern analyses in 
QCD based on the OPE such as e.g. the interesting work of Ref. 

The approximation which consists of substituting II(g 2 ) — > U OPE (q 2 ), with s large 
enough, has been given the name duality. Duality refers to the fact that, if the ap- 
proximation were exact and the substitution H(q 2 ) — > H.ope(q 2 ) carried no error, one 
could say that the integral of the experimental spectral function Im H(t) in Eq. (0) is 
dual to the OPE. The term duality violation, consequently, refers to any contribution 
missed by this substitution. 

The question is whether this approximation works for large sq, and if so, how large 
so should be for duality to apply within a precision of, let us say, 10%? In order to 
address this question, it is essential to be able to assess the validity of the approxi- 
mation used, lest the extraction of OPE parameters not be afflicted by uncontrolled 
systematic errors. Therefore, studying the contribution missed by the replacement 
II(g 2 ) — > XloPE^q 2 ) in Eq. (JTJ) is a necessary task in the present era of precision deter- 
mination of QCD parameters. 

That the OPE can at best only be an asymptotic expansion in the Euclidean region 
can be most neatly seen in the case of a two-point Green's function made of a pair 
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Figure 1: This figure shows the contour in the complex plane leading to Eq. f7)) . The crosses 
on the real axis denote the position of the resonance poles in the N c — ► oo limit, or the cut 
arising at finite N c . 

of quark bilinears, such as the vacuum polarization, in the large- N c limit [6]. This is 
because in this limit all meson resonances appear as intermediate physical states with a 
vanishing width and, consequently, as poles on the positive real axis in the q 2 complex 
plane. In contrast, the OPE is a series expansion in inverse powers of q 2 (up to loga- 
rithms) with no sign of a pole anywhere but at the origin, which means that it cannot 
converge on the positive real axis. The OPE can thus at best be asymptotic, 1 and even 
then only provided the positive real axis, i.e. the Minkowski region, is excluded. The 
region of validity of the expansion, therefore, is expected to be some sort of angular 
sector around the origin not including the positive real axis. Much less clear are more 
detailed properties such as, for instance, what the shape of this angular sector may be. 

Presently, the theory of duality violations in QCD is almost non-existing and it is 
important to develop it. In the absence of such a theory, we see no other choice but 
to resort to models. Somewhat surprisingly, not much work has been devoted to this 
issue. To our knowledge, Ref. [7] was the first to point out the importance of duality 
violations. That the OPE may also miss other important contributions has recently 
been suggested in Ref. jS]. 

In this work we will study a model of duality violations. We hasten to emphasize 
that our model is not QCD. In spite of this, we believe that it shares enough properties 
with QCD to be an interesting theoretical laboratory for studying issues related to 
duality violations. After all, our goal here is not a detailed numerical study of duality 
violations in QCD, but the much more modest one of understanding what the main 
features are, and how these violations may manifest themselves. This is not just an 
exercise of purely academic interest, however, because the lessons learned from the 
model raise questions which must also be addressed in the case of QCD. In addition, 
one of the results of our analysis will be a possible method for unravelling duality 
violations in QCD itself. 

1 A convergent expansion always has a circular region of convergence characterized by a non-zero 
radius. 
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The model starts from the Uyv ~^-aa vacuum polarization Green's function, in the 
chiral limit, given in terms of a spectrum consisting of two infinite towers of equally- 
spaced meson resonances, together with the rho meson and the massless pion, which 
are treated separately. Similar models can be found, e.g., in Refs. [71 1§1 ITTH ITU ITS]. 

The fact that the spectrum is infinite is crucial and underlies all the interesting 
properties to be discussed. This infiniteness of the spectrum is deeply rooted in the 
1/N C expansion of QCD. In fact, as a first step, we will take the large- iV c limit and 
consider a spectrum of resonances with vanishing width. This has the advantage that 
all quantities can be analytically calculated. As a second step, following Ref. [7j, we 
will give the resonances a width as a subleading 1/N C effect, while making sure that 
known analytical properties of the Green's function are not spoiled. We will thus be 
able to study a rather realistic model, allowing us to draw some conclusions on the 
issue of duality violations in the real world. 

It is obvious that, in order for Eq. (JIJ to be useful, an assumption has to be made. 
If the asymptotic regime starts at a value of q 2 which is higher than the sq employed 
in this equation, the use of the approximate expansion Hope{q 2 ) will completely miss 
all the physics of the full function IT(g 2 ), resulting in so large an error that, in fact, the 
approximation will become totally invalid ^3]. in the case of our model the asymptotic 
regime can only start after the contribution of the two towers of resonances has been 
taken into account, which means that the so scale has to be larger than the two lowest- 
lying states in these towers. The model might overemphasize the effects of duality 
violations, but, by the same token, it provides a means to identify the causes of these 
effects. 

Past experience accumulated over the years suggests that the case of real QCD may 
be less drastic than our model as far as the energy scale at which the OPE starts being 
a useful approximation is concerned. In other words, previous analyses of QCD which 
rely on the use of the OPE do not necessarily have to have systematic errors of the 
same size as those encountered in our model. But the model can be used to compare 
the different methods of analysis to see if there is one which is more reliable than the 
others. 

In real life, the values employed for so are in the interval ~ 1.5 — 3.0 GeV 2 and the 
lore is that this is expected to be high enough for the OPE to set in. The agreement 
within present errors in the determination of a s from tau decays [5] with those from 
other determinations at much larger energy scales [T3] can be taken as a confirmation 
of this. However, a much more accurate extraction in a new generation of precision 
experiments might show discrepancies between the different methods of analysis. In 
fact, there is already a glimpse of trouble in the determination of quark masses [To] , 
which becomes a more serious discrepancy in the terms which appear at higher order 
in the OPE, such as condensates |T0J . 2 We expect the lessons on duality violations we 
learn from our model to be of help in all these circumstances. 

2 For a nice summary of results, see Table 1 in Ref. ^3- Since certain electroweak matrix elements 
are related to condensates 0, one can also infer the latter from lattice determinations of the former. 
See Table 8 in Ref. [TBI - 
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2 Duality violations in the large- N c limit 

We will concentrate on the two-point functions of vector and axial- vector currents 

n^(g) = i I d'xe^(J^ A (x)Jt V > A (0)) = (q,q u - g, u q 2 ) U VA (q 2 ) , (2) 

where, for definiteness, Jy(x) = d(x) r y ,M u(x) and J A (x) = d(x) r y^ lr y 5 u(x). Throughout 
the paper, we will stay in the chiral limit for simplicity. Both functions, Uv,A{q 2 ) in 
Eq. (J2J) satisfy the dispersion relation 

^v,a(Q 2 ) = I 7 ^— ~ Imn^(t) , (3) 

Jo t - q 2 - It 7T 

up to one subtraction. 

According to the discussion in the introduction, and following Refs. [7J |5] , we will 
assume the following spectra: 

1 x 

-Imn v (i) = 2F 2 p 5{t-M 2 p )+2Y^F 2 {n)5{t-M 2 {n)) , 

n=0 

CO 

-lmU A (t) = 2F 2 5{t)+2Y j F 2 A (n)5{t-M 2 A (n)) . (4) 



7T 

n=0 



Here _F p is the electromagnetic decay constant of the p, and M p is its mass. Fy^{n) are 
the electromagnetic decay constants of resonance n in the vector and axial channels, 
while My t A{n) are the corresponding masses. Fq is the pion decay constant in the chiral 
limit. The dependence on the resonance index n is taken as follows: 

F vA n ) = F 2 = constant , M VA (n) = m 2 Y A + n A 2 , (5) 

which is compatible with known properties of the large- N c limit of QCD jH] , as well as 
properties of the associated Regge theory |T§] . 



The combination 



thus reads 



n(g 2 ) = ^(n v (g 2 )-n A (g 2 )) (6) 



Wo 2 ) = f± + l£ V - + V - (7) 

q^-qi + M* Z.- g2 + M 2 H +Z._ g2 + M 2 (n) • U 

Note that the upper end of the integration region in Eq. (j^J) requires the introduction 
of a cutoff, Aco, to render the integrals well defined. This cutoff must satisfy chiral 
symmetry in order not to introduce any explicit breaking through the regulator. This 
is accomplished by demanding that A^ ~ Ay A 2 ~ A^A 2 as N\r,A ~~ ¥ 00 • Of course, 
physical results are independent of the precise details as to how this limit is taken, in 
agreement with very general properties of field theory [TUj . 
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One can see this explicitly in Eq. (JJJ). The infinite sums are by itself divergent and 
have to be regularized. In order to obtain a convergent result, one has to choose the 
same regulator for both axial and vector towers, i.e. Ny = Na, in this case. Infinities 
then cancel and the two-point function can be expressed in terms of the Digamma 
function with the use of the identity 

N , N 



E^-E 1 --7*, (8) 



lim 

n^qo I z — ' z + n n I v ' z 

. n=l n=l 



where 7^ is the Euler-Mascheroni constant and ip{z) is the Digamma function, defined 

as 

= / *=5;^ W . (9) 



Thus, it follows that 



_p2 _p2 T / ,;l ,,,2 \ / ,,2 ,,,2 



n ^ ^ g 2 + -g 2 + M 2 + A 2 



— q + m A \ I —q + m v 



a 2 y V A2 

Furthermore, we will impose the following conditions on our model: 

2 JV„ F 2 



(10) 



316vr 2 A 2 ' 

so that the coefficient of the parton-model logarithm is given by a free-quark loop; 



(11) 



2 1\ „, / m 2 A 1 



' ' ['-j^-i) > ^ = ^ 2 ^--) (12) 



2 U 2 / . / "V m V , 1 \ r 2 a 2 / m A m A 1 



and 

-2F^ + F^^-^ + -J=W^-^ + ^ , (13) 

so that, as in the case of real QCD, there are effectively no dimension-two or dimension- 
four operators in the OPE of IT(g 2 ) (see below). Note that Eqs. (JT2J) and (fT3*j) would 
not be consistent if F p = 0. Since we want our model to be semi-realistic, this is the 
reason why we separated the rho meson from the vector tower in Eq. (J3J. This is also 
natural from the point of view of Regge theory, as the rho-meson mass is low enough 
to be below the onset of Regge trajectories (presumably at a scale of 0{1 GeV)). The 
di is heavier and thus we include it in the axial tower. 3 At any rate, we will see that 
what really matters is the analytic properties of the OPE and these depend on the 
tower as a whole, rather than on the properties of any individual resonance. 

Knowledge of the spectrum (0J) allows us to calculate the full Green's function 
II(g 2 ) explicitly. Defining the spectral function p(t) = (1/tt) Im II(t) one has that 
(with Q 2 = -g 2 ) 

POO J 1 f~i 

= / f^-*) A £ p, (14) 

^ u ^ fc=l,2,3,... ^ 



3 We could also have separated the a%, but that would introduce unnecessary free parameters 
without any crucial change in the analytic properties of the OPE. 
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where the 1/Q 2 expansion in the second equation above is akin to the Operator Product 
Expansion in QCD. Note that a naive expansion in powers of t/Q 2 under the integral 
sign is not an allowed mathematical operation. Consequently, in order to obtain the 
correct 1/Q 2 expansion due care must be exercised. The correct OPE can be derived 
by using (the derivative of) Eq. Q and the expansion 



oo 1 

e z - 1 ^ y n\ 

n=0 



(15) 



wherein B n (x) are the Bernoulli polynomials. In this way one arrives at the following 
expressions for the OPE coefficients C^k- 



a 



21,- 



-F 2 S k>1 + 



-1 



F 2 M 2k-2 



2 a 2fc-2 



k' 



-FA 



m 



A 2 



In particular, the first few coefficients read 4 



C 2 
C 4 

c 6 



—F 2 M 2 + -F 2 A 2 



Bo 



+F 2 M 4 p - -F 2 A 4 ■> B, 



-F 2 p M % p + -F 2 A 6 



5, 



and the first few Bernoulli polynomials are 



B (x) = l , B 1 (x)=x- 



1 



B s (x) 



n ^9 1 

x x H — x 

2 2 



2 1 
B 4 (x) 



Bi 



my 

A 2 

m| 

A 2 

m v 

A 2 

m v 

A 2 



B 2 (x) 



B k 



m , 



V A 2 



(16) 



Bo 



-B. 



B d 



m A 

A 2 

m A 

A 2 



m , 



\ A 2 



(17) 



x — x + 



X 



2x 3 + x 2 



1 

6 
1 
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As already stated, the OPE cannot be obtained by naively expanding at large Q 2 
under the integral sign in Eq. (jH| . We emphasize that this feature is not specific to 
our model and is, in fact, also true in QCD. Note that, if a naive 1/Q 2 expansion 
were valid in QCD, there could be no logQ 2 dependence in C<ik and all anomalous 
dimensions would have to vanish, which is obviously not the case. Our model does 
share with QCD the property of a non-trivial expansion in 1/Q 2 but, in the large- N c 
limit, it still misses all these log's. With the introduction of finite widths in the next 
section, this drawback will be somewhat ameliorated. 

It thus follows from the previous discussion that the moments 



MJs ) 



AO 



dt t n p(t) 



(19) 



4 In fact, C24 vanish identically because of the conditions 
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are not the coefficients of the OPE, C 2 k- In order to obtain the actual relation be- 
tween M n (s ) and C 2 )., it is convenient to first calculate the Laplace transform of the 
dispersion relation (|14j) 



poo poo 

n(-g 2 ) = dr e- rQ2 p(r) with p(r) = dt e~ tT p(t) . 
Jo Jo 



(20) 



Plugging in the 1/Q 2 expansion of Eq. (j!4|) and identifying Laplace transforms on both 
sides, one obtains 

oo n 

^2k 



E 



.k-l 



dt e~ tT p(t) 



and therefore, finally, 



a 



2k 



lim { 



fc-i 



(21) 



(22) 



As an exercise one can use the expression (|22|) to recover our results in Eq. ()17|) . It 
is the slow fall-off with t of p(t) which is responsible for all the complications and, in 
particular, the duality violations. If p(t) were falling off faster than any power at large 
t, one could take the r — > inside the integral, and the naive 1/Q 2 expansion of the 
dispersion integral in Eq. (|14|) would be valid. 

In actual fact, the calculation of the moments (fTU|) yields 



Mo(so) 

M x (s ) 
M 2 (s ) 



C 2 -F 2 
-C 4 - F 2 



C 6 - 

— -F 2 A 4 
3 

-C 8 -F 
-F 2 A 4 



B x {x v ) - Bi(x A ) 

B^Xy) - B x {x A ) 
Bi(x v ) - Bi(x A ) 
B 3 (x v ) - B 3 (x A ) 
Bi(x v ) - B x {x A 
B 3 (x v ) - B 3 (x A ) 



so + ^F 2 A 2 



B 2 (x v ) - B 2 (x A ) 



F 2 A 2 



B 2 (x v ) - B 2 (x A ) 



so 



sl + -F 2 A 2 [B 2 (x v )-B 2 (x A ) 



so + -F 2 A 6 



B A {x v ) - B 4 (x A ) 



(23) 



where the C 2 k are given in Eqs. (|T7j) . and we have defined < xy, A < 1 as the fractional 
parts of (s — m^)/A 2 , respectively. The above expressions are valid when s > m v , 
i.e., when both the vector and axial towers are included. In these expressions, the 
points xy >A = 0, 1 should be excluded because they correspond to the singularities of 
the delta functions in the spectrum, cf. Eq. (0J). We have kept the coefficients C 2A 
for illustrative purposes, even though they actually vanish because of the constraints 
(jlHll3j) . Perhaps several comments to explain the pattern in the result of Eq. (J23)) are 
in order. 

First, these expressions clearly show that there is no So at which the moments 
equal the OPE coefficients, C 2 k, and consequently, in the world of N c — > oo described 
by our model, there is no such thing as a true "duality point." There are always 
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Figure 2: From left to right: Mq, M\ (first row), M 2 and M3 (second row), defined in 
Eq. iiy\) . as functions of the energy scale, sq (in GeV 2 ). 



duality violations in the form of a polynomial in the variables xy,A- As can be seen in 
Fig. these violations are step-function "oscillations" around the corresponding OPE 
coefficient. Note that the OPE coefficients Ce,s are numerically negligible (using the 
parameters in Eq. (|24|) below) compared to the oscillations shown in Fig. |2] for the 
corresponding moments M 2 $. 

Looking at Eqs. (}2"3*j) . one sees that there is some recursiveness, i.e. higher moments 
depend on the duality violation already existing in lower moments. Also, higher mo- 
ments have higher powers of Sq modulating their duality violations, so the oscillations 
grow larger for higher moments. This means that the larger the Sq at which one is 
computing the moments the worse the duality violation. All this behavior is 

clearly shown in Fig. where we have taken as the values for our parameters the set 5 

F = 85.8 MeV , F p = 133.884 MeV , F = 143.758 MeV , (24) 
M p = 0.767 GeV, m A = 1.182 GeV, m v = 1.49371 GeV , A = 1.2774 GeV . 

Second, as we have said, the end points xy,A = 0, 1 are to be excluded from the 
result shown in Eq. (|23|). They correspond to the singularities of the delta functions in 
Eq. (j3J), i.e. to the position of the meson masses, and cause the appearance of step-like 
discontinuities at the same location in all the moments M n (so). At the jump produced 
by, e.g., the kick of a vector resonance after the drop produced by an axial- vector one, 
Fig. |21 shows how one crosses the value of the corresponding coefficient of the OPE. 
However, since the slope at these points is infinite, there is no way to predict the C<ik 
coefficient from the value of the corresponding moment. 

5 This choice is dictated by the fact that, when a nonvanishing width will be included in the next 
section, the spectral function resembles qualitatively the experimental data. See Fig. below. This 
choice of parameters also satisfies Eqs. 
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One might think that the introduction of a finite width could drastically improve the 
present situation at large N c . However, this is not really true in the model. Although 
introducing a finite width certainly changes the slope to make it finite, it also moves 
the location of the duality points from one moment to the next. As we will discuss in 
more detail in the next section, the net result is that it is still very difficult to make 
accurate predictions. 

The best way to try to understand the relationship between the moments and 
the coefficients of the OPE is through the use of Cauchy's theorem, Eq. ([[J. As 
already discussed in the introduction, the replacement n(g 2 ) — > Hqpe^ 2 ) would be 
valid if the OPE were a convergent expansion with a radius of convergence including 
the circumference \q 2 \ = Sq. However, in QCD the OPE is not a convergent expansion 
anywhere on the Re q 2 > axis. Since this crucial property is shared by our model, we 
can now use the model to study how Eq. JI} should be modified when the expansion 
^-ope(q 2 ) is used instead of the exact function, H(q 2 ). 

To this end, let us define V^(sq) from the equation 



[ S0 dt t n pit) = / dq 2 q 2n U OPE (q 2 ) + P [ 

J0 2m J\q 2 \=s a 



[so) 



(25) 



Obviously, the term V^(sq) measures the amount of duality violation and we would 
like to obtain a more explicit expression for it in our model. 

In order to do this, note that the OPE in our model is related to the large \z\ 
expansion of the Digamma function, 



j oo 

ip(z) ~ logz - — - 



B 



2n 



71=1 



In z 



2n 



arg z < TT 



(26) 



which is not a valid expansion on the negative real axis in the complex z plane. We 
may, nevertheless, obtain a valid expansion in this region if we use the "reflection" 
property of the Digamma function 



1p(z) = lj){ — z) — 7T COtijtz) 



(27) 



which is valid for any complex number z. Using Eq. (|27jl . the appropriate expansions 
for large values of \q 2 \ in the whole q 2 plane then read 



, .TWO? 2 ) + O (e- 2 ^ 2 l/ A2 ) ,Reg 2 <0; 
\UoP E (q 2 ) + A 00 (g 2 ) + O (e- 2 ^ 2 l/ A2 J , Re q 2 > , 



where 



Aoo(g 2 



ttF 2 
A 2 



cot 



7T- 



-q 2 + m^ 1 
~A~ 2 



— cot 



TT 



-q 2 + m 2 ^ 
A~ 2 



(29) 



and the subscript "oo" refers to the N c — > oo limit. Let us emphasize again that the 
appearance of Aoo(g 2 ) in Eq. gives rise to a violation of duality which does not 
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Figure 3: Contour used for the definition of the functions D^ cin (sq) and T>^ nst (so) in 

Ens, imm . 

go away as \q 2 \ — > oo (in particular it contains simple poles on the positive real axis 
in the q 2 plane). This duality violation comes on top of the expected exponentially 
suppressed contributions, which we have exhibited explicitly in Eq. ()28j) . and which 
originate from the fact that the OPE, in those regions where it is valid, is only an 
asymptotic expansion. 6 

Consequently, up to an exponentially small contribution, the duality violation T> 
term can be written as an integral over the semicircle \q 2 \ = so with Re q 2 > (taken 
counter-clockwise) of the function A^g 2 ), 

PW(So) = ~2^ / ^ ^ A °° (g2) + ° ( e ~ WA2 ) • (30) 
\<?\ = so 
Re q 2 > 

Furthermore, Eq. (J3U|) can be rewritten in the form of a spectral integral in the 
following way. The integral over the contour in Fig. El vanishes, since it encloses no 
singularity. It follows that the duality violation function T)^(so) in Eq. (|3U|) can also 
be expressed, up to exponentially small terms Q(e~ 2ns °/ A ), as 

^ ln \so)^V [ :l iU iso)+V^ st (s ), (31) 

with 

r i r so 1 

V [ :Lu.(so) = / dtr-lmA^t + ie) (32) 
V [ ±t.(s ) = ~[J tE + £°] dq 2 q 2n A^q 2 ) . (33) 

In these expressions, the integral fl32(l is taken over the positive real axis, in the complex 
q 2 plane, whereas the integrals ()33)1 are taken over the imaginary axis (with e — > + ). 

6 See the Appendix. 
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We see that expression (p£2|) is of the form of a spectral integral, just like the left-hand 
side of Eq. (J25j) . In fact, for N c — > oo, Im A QO (g 2 + ie) is also a sum of delta functions 
with poles precisely located according to the spectrum given by the towers of Eqs. (J5J). 
This follows from the expansion for the cotangent, 

1 °° 1 

7rcot(^) = - + 2^— -. (34) 

n=l 

Physically, the function T>^J cill (sq) is responsible for balancing all the oscillations 
which are present in the spectral integral on the left-hand side of Eq. (J2*5|) but which 
are missed by the contour integral of the OPE on the right-hand side. The function 
^coLi.( s o) ; 011 the other hand, does not have any oscillating component, and is domi- 
nated by the region of small q 2 , decaying exponentially fast in sq to a constant. One 
has 

2>£L(*o) = Cg+ 0(e- 2 -°/ A2 ) , (35) 
with, in our N c = 00 model, 7 

Using this result for n = 0,1,2,3 in Eq. flU together with Eqs. flHl, (JED and (jSHJ) 
one may recover the results of Eq. (J2*H|) for the moments Mo,!^- Again, the total 
contribution to T>^(so) in Eq. (f23J) does not vanish, even in the so — > 00 limit. 

Since duality violations are caused by the lack of convergence near the positive 
real axis in the q 2 plane, a possible strategy to eliminate them might be to use suitable 
polynomials Pit) in Eq. ()25|) which suppress the contribution precisely from this region 
in q 2 by having a zero at q 2 = Sq. Such polynomials have been referred to as "pinched 
weights" in Ref. [H]. 

This is confirmed by our large-Ac model, and it is related to the recursiveness 
present in the different moments in Eq. (|23|) . For instance, a pinched weight of the 
form (1 — t/so) automatically cancels the duality violating oscillation proportional 
to E>x{xy) — Bi{x a) between the moments M and M\. The cancelation, however, 
is not complete and this pinched weight leaves behind an oscillation proportional to 
[B 2 (xy)— B 2 (xa )]/sq which, although suppressed at high s , does not vanish completely. 

One can actually exploit this feature a little further and, by inspection of Eq. (J23J), 
one easily sees that P 3 = (1 — t/s ) 3 is the pinched weight which best suppresses the 
duality violating oscillations among the moments M^i 2,3 since it only leaves behind a 
term proportional to \B^{xy) — B^(x a)]/ s^. This is illustrated in Fig. EJ where it is 
shown how much of a reduction in the oscillations of the left-hand side of Eq. (J25|) is 
achieved by the pinched weight P3 in comparison with other pinched weights previously 
proposed in the literature i-e. P\ = (1— 3t/so)(l— t/so) 2 and P 2 = (t/so)(l—t/so) 2 . 

Figure H] also shows that the OPE curve more or less interpolates through the 
oscillations in the spectral integral. As one can see, there are points at which the 

2 2 

7 In order to obtain this result, note that in our model 1 < ^ < 2, whereas ^ < 1, cf. Eq. J23J). 
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Figure 4: Integrals of the pinched weights, Jq° dt -Pi,2,3(£)p(£) defined in the text (solid lines, 
from left to right and top to bottom) and the corresponding OPE prediction (dashed lines) as 
a function of sq (in GeV 2 ). These integrals are expressed in units of 10 -2 GeV 2 . 



oscillation happens to cross the OPE curve so that, at these points, all duality violations 
effectively vanish. Based on this, one could imagine guessing the OPE coefficients by 
some sort of averaging over one full oscillation. Regretfully, this behavior sensitively 
depends on the fact that we are considering the N c — > oo limit. The inclusion of a 
finite width essentially wipes out any oscillating signal and one is left without any clue 
as to where to take the average. 

In the next section we will discuss in detail the inclusion of a finite width as well 
as those features which survive the passage to this more realistic setting. 



3 Duality violations at finite N c 

Since the issue of duality violations hinges strongly on the analytic properties of the 
OPE and the corresponding Green's functions, it is clear that the introduction of a finite 
width has to be done with some care in order not to spoil these analytic properties. 

Moreover, in QCD finite widths appear as a consequence of subleading effects in 
the 1/N C expansion. It is clear that if one knew how to properly include all these 
1/N C effects, one would incorporate just the right 1/N C corrections to give resonances 
a finite width, while at the same time preserving the relevant analytic properties of 
Il(g 2 ). Among those properties, the most important one for us will be the existence 
of a cut in the complex plane only for Re q 2 > 0. A simple Breit-Wigner function, for 
instance, is not good enough. 

However, the solution of QCD in the large- N c limit is unknown, let alone the 1/N C 
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corrections, so it is clear that one has to model these effects in some simple way in 
order to make progress. This is what has been done in Ref. [7]. The hope is that 
the conclusions drawn from such a model will not depend in any crucial way on the 
particular details but rather on the generic properties of the model. 

With this philosophy in mind we follow Ref. [7f and introduce a subleading l/N c 
width, writing the vacuum polarization H(q 2 ) as 



irNj z z + M* ^z + M\{n) ^z + M 2 (n) 

r n=0 n=0 v 



(37) 



where now the variable z denotes the combination 



and My^(n) are given in Eq. ( |5|) . As in the zero-width case, the infinite sums can be 
expressed in terms of the Digamma function, and one obtains 



1 f F 2 F 2 F 2 



. z + m 2 A \ ( z + m 2 v 



(39) 



Obviously, in the iV c — > oo limit we recover the results of the previous section. The 
conditions of Eqs. (fTT|) , (|T^jl and (fT3~j) impose now the absence of any log z,l/z and 1 / z 2 
terms in the large Q 2 = —q 2 > expansion of the Green's function H(q 2 ) in Eq. (|39p. 
The choice of values for the parameters in Eq. (|24|) plus the choice 



a = 0.72 (40) 

for the new parameter a produces the spectral function which is plotted in Fig. |S| 
As can be seen in this figure, the model is not able to reproduce the experimental 
data in detail, but the qualitative agreement is good enough for our purposes. We 
remind the reader that our model is not meant to extract physical parameters, but 
only as a laboratory for testing different methods which are currently in use for the 
determination of the coefficients of the OPE from the data. 

In the limit N c — > oo one recovers for Im n(g 2 ) the result with Dirac deltas of the 
previous section, Eq. (jlj. However, at small but finite a/N c , the dependence of every 
resonance propagator on q 2 through the new variable z produces a complex pole with 
an imaginary part, and thus a width Tv,a(ji) given by 

TvAn) = ^-MyAn) {} + ° {jfj) > ( 41 ) 

and, mutatis mutandis, for the p as well. This behavior for T is the one expected on 
general grounds jTj and, as one can see, forces T VA (n) to vanish as N c — > oo, although 
it should be noted that the width also grows with the resonance number n. 
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Figure 5: Comparison of the data from ALEPH and OPAL (overlaid on a single plot and 
including the kinematical factors of r decay) with the model, ^Im H(q 2 ) from Eq. \3y\) (solid 
blue line), as a function of q 2 (in GeV 2 ). As it is the case with the experimental data, the 
pion contribution, Fq/z in Eq. i,^y\) has been subtracted away. 



Given that the change produced by the widths amounts to the replacement 
the expansion of the Green's function n(g 2 ) at large values of the momentum \q 2 \ is 
similar to that in Eq. (}2T 



n(g 2 ) 



except that now 



n 



OPE 



O e 



-27r|g 2 /A 2 | C 



Jl PE(q 2 ) + A(q 2 ) + 0(e- 



■2w\q 2 /A 



2|C 



Re q 2 < , 
Re q 2 > , 



(42) 



A(g 2 ) 




7T 



A 2 



+ 7T 



A 2 



— cot 



7T 



A 2 



+ 7T 



A 2 



(43) 



with ( = 1 — -j^. Similarly, the OPE is given by an expansion as in Eq. (fPfj) but the 
coefficients (fTBj) are shifted by a log Q 2 -dependent term. For small a/N c , the Wilson 
coefficients read 



C*2fc(Q 



2\N C =3 



c»r*> 1 1 + 



a ka Q 2 ( a 2 \ 



(44) 



where the C 2fc c 00 are those given in Eq. (|16p. We will find it convenient to explicitly 
separate the log dependence and decompose these coefficients as [20] 



C 2k (Q 2 ) 



,2\7V C =3 



Q 2 ( a 2 
a 2k + b 2k + 



(45) 



The b 2k would model the existence of anomalous dimensions in real QCD were it 
not for the fact that the latter, unlike the b 2kl do not vanish in the large- N c limit. 
We do not consider this to be a serious drawback of the model, however, and the 
1/N C contributions in Eq. (j44j) still mimic the a s (m T ) contributions to the anomalous 
dimensions which exist in real QCD. Even though the 1/N C corrections in Eq. (j4^j) 
grow larger with increasing dimension of the operator, at least for the first terms there 
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is a sense in which the N c — * oo limit is close to the world at N c = 3 in the OPE. 
Plugging in our set of numbers (|241 14U|) . the OPE coefficients are 

a-2 — i>2 = = a 4 = 6 4 
a 6 = -2.8 x 10~ 3 GeV 6 , b % = -5.9 x l(T 4 GeV 6 , 
a 8 = +1.8 x l(T 3 GeV 8 , b 8 = +5.1 x 10~ 4 GeV 8 . (46) 

The small 1/N C correction to the Green's function II(g 2 ) in the Euclidean (cf. 
Eq. (|44|) ) is in complete contradistinction to what happens to the spectral function 
Im n(t) in the Minkowski region, t > 0. As emphasized in Ref. [7j, the reason for 
this is that the limits N c — > 00 and t — > 00 do not commute in this region. Defining 
q 2 = t > 0, one obtains to leading order in a/N c for the variable z of Eq. (J38|) in this 
region, 

a , ( t \ .a _ / a 2 



-- ( ( 1 ^ ,og UJ + ^ + °llv?JJ ' (47) 

Taking imaginary parts of the reflection property for the ip function in Eq. ()27jl . at 
large values of the momentum t > one obtains that 



where 



and 



Im II(t) = Im n Pi?(t) + Im A(t) , (4J 



Im A(t) = — c sr sm ^ ^ Kj sm ^ ^ (^1 + ^- 

+ O fe" — ^) . (50) 

The imaginary part in Eq. (J4T?|) can be gotten from the imaginary part of the logarithm 
in the OPE expansion in Eq. (|44j) and is also familiar from the QCD case, except for 
the 1/N C suppression. Analogously, the cancelation of the t" 1 and t~ 2 terms is due to 
the absence of the corresponding 1/Q 2 and 1/Q 4 terms in the OPE in the Euclidean 
region. 

The contribution shown in Eq. (|5Uj) is completely missed by the OPE and stems 
from the duality- violating function A(g 2 ) in Eq. (|43|) . We call the reader's attention 
to the N c factor in the denominator of the exponent in Eq. (|50|). The expression in 
this equation is obtained in the large t limit when A^ c is kept finite. That is to say, for 
finite N c , the isolated poles of Aoo(g 2 ) on the real axis (Fig. EJ) become a cut. If, on the 
contrary, the limit N c —>■ 00 is taken first, then the expression (|5U|) is not valid and one 
obtains the infinite sum of delta functions of the previous section instead. It follows 
that the two limits N c —>■ 00 and t — » 00 do not commute in the Minkowski region. We 
expect this effect to be quite generic, and we thus believe the exponentially damped 
oscillation shown in Eq. (|50j) to be more general than our particular model [7]. 
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In general, although Eqs. ()25I30|) are still valid, the behavior of the duality- violating 
function A(g 2 ) changes drastically with respect to its counterpart A 00 (g 2 ), showing an 
exponential fall-off at large \q 2 \ on the half plane Re q 2 > 0. One has now that 



2tt \ ^hr\ sin 



We see in particular that the limit iV c — ► oo (i.e. £ — *■ 1) makes the values (p = 
0, 27r exceptional: when N c = oo, the exponential fall-off on the real axis completely 
disappears. This is of course consistent with Eq. ([5(J|h 

The exponential fall-off exhibited in Eq. (J51|) implies now that the function V^(sq) 
vanishes exponentially at large sq as 

D [n] (s ) ~ O (e~^^) . (52) 

Again we note the singular limit N c — > oo in this expression: in the previous section 
we saw that T>^ n \So) —> "undamped oscillations," if the limit N c — > oo is taken before 
So is taken to infinity. 
For finite N c , one has 

= ilW+CW ^ O(e-^B), (53) 
2>£Lh.(«o) = rdtt n UmA(t + ie) Ae ~C [n] + O (e- 2 ^) , (54) 

J 

Z>£L.(«o) = CW + (e~ 2 ^) . (55) 



We are using the hierarchy 

e" — > e ~ 2 ^ j (56) 

which is very well satisfied for large values of Sq. Note the more stringent bound on 
the corrections to Eq. ([51)]). which results from the fact that T>^} nst (s ) is defined as an 
integral over A(q 2 ) for values of q 2 far away from the positive real axis. 8 
Combining Eqs. (j53H55j) in the sq — > oo limit, it follows that 



poo 1 

CN = _ / dtt n -Im A(t + is) , (57) 
Jo 7r 

and, inserting this result into Eq. (|53|). we obtain 

dt t n - Im A(t + ie) + O (e-^) (58) 



as the final expression in the finite width case, for the duality violations in Eq. 
Since Im A(t + is) ~ C^e" 1 ^^ 2 ") (cf. Eq. (|50|)). the integral in Eq. (JHEJ) is indeed of 

27TSQ 

order e ) in agreement with Eq. (|53j). 
8 See the Appendix for the derivation of Eq. (pSJ. 
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Figure 6: From left to right: Mq, Mi (solid curves, first row), M2 and M3 (solid curves, 
second row), defined in Eq. MfJ\) . together with the OPE prediction (dashed curves) as given 
by the corresponding terms in squared brackets in Eqs. I5yii6ty) . as functions of the energy 
scale, sq (in GeV 2 ). 



Returning to Eq. 1)25)1. we obtain for the OPE coefficients the following set of 
equations: 

h b 8 



SO 

dt p(t) + 


SO 

dt t p(t) + 



2sq 3sq 



+ 



+ 



dt t 2 p(t) - 



a 6 + b e log — H h 

A 2 s 



dt t 3 p(t) + 



■ 7 1 s o , bio 

a 8 + 6 8 log — - b 6 s H 

A 2 s 



^ [1] (^o) 



(59) 
(60) 
(61) 
(62) 



so 



dt t 7 p(t) + 



a M + 6ie log ^ 



^8 g Q _|_ ^18 

4 "' s 



rt 7 \s ) , (63) 



where the right-hand side is given by Eq. ([58)1 . The results for the lowest moments 
()19|) and the corresponding OPE (the terms explicitly shown in squared brackets in 
Eqs. (J2H2D) are plotted in Fig. El 

We draw the following conclusions from Eqs. (|59H63)) : First, it is clear that the 
b coefficients of the OPE also contribute to these equations. The pattern is simple 
and governed by dimensional analysis: a higher order b contributes to the equation 
determining a lower order a accompanied by inverse powers of the scale Sq. On the 
other hand, a lower order b contributes to the equation determining a higher order a 
accompanied by positive powers of s . 
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Even though in real QCD the b coefficients are order a s effects with respect to the 
a's, the fact that they may appear multiplied by positive powers of s tells us that, in 
general, it is safer not to neglect them from the start, unlike what is done at present 
in common practice. A possible exception is the term , which is the only one which 
has been calculated [2T] and which seems to be small enough to be safely disregarded. 
Since this term is the only one accompanied by positive powers of s in Eqs. (|59H62jl . 
it might not be too unreasonable to neglect the b terms in these equations in a first 
approximation. However, for higher moments, this practice may be dangerous. Notice 
that an Sq term amounts to a factor of almost 100 GeV 8 for sq = w? T = 3.15 GeV 2 so 
that, for instance, the term 6 8 Sg in the last Eq. might not be negligible even if b 8 
is small. 

Second, and most important, the right-hand sides of Eqs. (J59H63)) do not vanish in 
general, except in the limit So — > oo. As it is obvious from these equations, the D 
terms may potentially pollute the extraction of the OPE coefficients. Consequently, 
in any precise determination of OPE coefficients it is unavoidable to take this duality 
violation into account. 

In order to further appreciate this point, let us close our eyes, neglect the duality 
violations T> altogether, and apply the different methods which have been employed in 
the literature so far. 

For instance, let us start with finite-energy sum rules. The method consists in 
determining a duality point Sq from the condition M .i(sq) ~ 0, which is then to be 
used in M 2j 3(,So) to predict the combinations appearing in Eqs. (j61|62|) . namely 

4(s*) = a 6 + 6 6 log^| + ^ + ... , 

A 8 (s*) ee a 8 + 6 8 log^|-6 6S * + ^ + ... (64) 

A S 

In our model the first duality point happens to sit at Sq = 1.472 GeV 2 , yielding the 
predictions 

A FESR = _ A q . 1Q -3 Ge y6 ; A FESR = +Q 3 . 1Q -3 Ge y8 ^ (65) 

which can be compared to the correct values 

45(1.472 GeV 2 ) = -2.4 • 10~ 3 GeV 6 , 4,(1.472 GeV 2 ) = +2.6 • lO^GeV 8 . (66) 

There is a second duality point around 2.4 GeV 2 fulfilling either M (sq) = or Mi(sg) = 
0, but not both. Choosing the one which satisfies Mi(sq) = 0, for instance, one gets 
s* = 2.363 GeV 2 and this yields 

A% ESR — —2.0- 10~ 3 GeV 6 , Af^ = -1.6-KT 3 GeV 8 , (67) 

to be compared to 

4(2.363 GeV 2 ) = -2.8 ■ 10" 3 GeV 6 , 4(2.363 GeV 2 ) = +3.4 ■ 10" 3 GeV 8 . (68) 
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Figure 7: Difference between the integrals of the pinched weights, Jq° dt Pi t 2,3(t)p(t) and the 
corresponding OPE contribution from Eqs. \59\6 6 A) (in units of 10 -3 GeV 2 ), as a function of 
sq (in GeV 2 ,), from left to right and top to bottom, respectively. 

As we see, in going from the first to the second duality point, A 6 is reduced by a 
factor 2 and A 8 changes sign. Interestingly enough, this is precisely the trend seen in 
the corresponding determinations in the literature [22 I2H] • 

From the comparison of both results and (|67|). it seems that there is no dramatic 
advantage to making the prediction using the second duality point rather than the first, 
even though at the second point Sq has a larger value. 9 Moreover, while going to a 
larger s is m principle better because the duality violations T>^(sq) become smaller, 
if this is done while keeping all the b coefficients to zero, i.e. taking the A's above as 
estimates for the a's, then the growth of the moments M2${sq) with sq — they are 
actually divergent because of Eq. (|49|) — is numerically absorbed by a wrong shift in 
06,8- As we see in Eqs. and (|57|). this may invalidate the advantage of choosing a 
larger s - Furthermore, in the real world the experimental error bars are much bigger 
at larger sq, the QCD value of is very small and that of 6§ is unknown. This makes 
it difficult to know exactly what explains the difference between the results obtained 
in the literature using finite-energy sum rules [22 123 El- 10 

Let us next consider pinched weights, which have been suggested [201 I2S| as a 
way to reduce the amount of duality violations. However, the use of pinched weights 
is no guarantee that one can do away with the duality violations altogether. For 
instance, we have followed Ref. [20] and fitted the OPE coefficients to the combination 
of moments obtained with the pinched weights P± = (1 — 3t/so)(l — t/so) 2 and P2 = 

9 The prediction for A§ at the second Sq is somewhat better than at the first, but this is no longer 
the case for A$. 

10 Reference [21] finds numerical agreement with Ref. [23 but, in fact, it does not have the second 
duality point the analysis of Ref. [231 is based on. 
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(t/s )(l — t/s ) 2 , neglecting duality violations. We have done this with 20 points in the 
window 1.5 GeV 2 < s < 3.5 GeV 2 . The upper end in this window was identified by 
noting that the tau mass roughly coincides with the mass of the p(1700) in QCD, and 
in our model this resonance happens to be at m 2 „ ~ 3.8 GeV 2 . The result from the fit 
can be interpreted as some sort of average over the window in so of the combinations 
A 6i8 in Eq. (|51jl. and yields 

A pinch = _ 3 g x 1Q -3 Ge y6 ^ A pinch = g g x 1Q _ 3 Ge y8 _ (69) 

As one can see, the error made is again rather large (particularly in the case of As) 
and, in fact, comparable to that using finite-energy sum rules in Eqs. (|531I57|) . We did 
not see any clear improvement in precision by changing details of the fit such as, e.g., 
the window of s Q values used. 

Since pinched weights suppress duality violations, one may try to design "good" 
pinched weights which have as small duality violations as possible. Assuming a general 
expression for the duality violations of the form (cf. Eq. (JSQ]) ) 

-ImA(t) = k e-^ sin (a + fit) , (70) 

7T 

for certain values of the parameters k, 7, a and (3, one could ask which is the pinched 
weight, involving moments not higher than M 3 , which minimizes the amount of duality 
violations. As before, the answer to this question is P3 = (1 — t/so) 3 . The reason is 
that, given the general form in Eq. (|7())1. one finds that P^s residual oscillation is 
damped by the factor So 3 e _7S °, while in the case of P\^ the damping factor is only 
So 2 e _7S °. This is a behavior reminiscent of the N c = 00 case in the previous section. 11 
Figure [7| shows the difference between the integral of the pinched weights Pi 2,3 an d 
their corresponding OPE contribution from Eqs. (j59H62j) as a function of Sq. If there 
were no duality violations these curves should be a flat zero. As this figure shows, P3 
has the smallest duality violations. 

In Ref. [17J a Laplace sum rules analysis was done. Neglecting the b coefficients one 
obtains for the a coefficients of the OPE the equations 

,L „L 1 _ / j. -tr „/.\ , Z / j. . -tr 



«6-aioT^ + - = — / dte- tT p(t) + - dt t e~ tT p(t) , 
11 T Jo T Jo 

t 12 f s ° 6 f s ° 

4 + af - + ... = -— / dte- tT p(t)-- dtte- tT p(t). (71) 

Requiring stability under variations in the Laplace variable r one can determine ag 8 . 12 
In Fig. |H1 we have plotted the true values for a 6>8 from Eq. together with the r 
dependence of ag , a 8 - neglecting af and higher. The upper limit in the Laplace 
integral has been chosen to be the second duality point from the previous finite-energy 
sum rules analysis of Eq. (JHZJ), Sq ~ 2.35 GeV 2 , but there is no qualitative change if the 
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This is a general result. For moments not higher than M n the optimal pinched weight is (1— t/so) 



-7S0 



and the residual oscillation is modulated by s n e 

12 Obviously, the true a coefficients (|46H are independent of r. 
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Figure 8: Plot of the coefficients a§ (solid curve, left) and Og (solid curve, right), together 
with the corresponding true values (dashed horizontal lines) from Eqs. (all in units 
of 10~ s ), as a function of the Laplace variable, r (in GeV~ 2 ). 



first duality point is used instead. As it can be clearly seen, no sign of stability in 
r is found. (We do not expect this to change qualitatively if a w , b 6 etc. are included.) 
This is unlike what is found in QCD and may be another sign of what we already 
said in the introduction concerning the model's tendency to maximize the violations 
of duality. 

Finally, we also used the MHA method [26 . In this case one truncates the spectrum 
to a finite number of resonances whose parameters are adjusted so that the low- and 
high-energy expansions of the II(g 2 ) function in the large-iV c limit are reproduced. In 
this way one constructs an interpolator for the function II(g 2 ) , called U.MHA(q 2 ) with 
which it is possible to calculate, e.g., the integrals which determine the relevant cou- 
plings governing semileptonic kaon decays. However, in the present case the situation 
is slightly different because one wishes to use ITm/m(<? 2 ) to predict coefficients of the 
high-energy expansion, and not as an interpolator for an integrand. 

The function Hmha((1 2 ) is constructed as 13 

H(q 2 ) MHA = % + ^JZ . — , (72) 

Q My — q 2 M\ — q 2 

with the parameters My, Ma, Fy, Fa to be determined from the conditions 6*2,4 = 
and the values of L 10 and the pion electromagnetic mass difference in the large-iV c 
limit, as taken from the model in the previous section, Eq. (JTUJ). 
They are given by the following expressions: 

LlQ = A F 2 M ± + Mv = _ (5 2 ± Q 5) . 1Q -3 (73) 
4 m 2 M\ K ' 

K - K = (^-,) Fl3^L log (SI ) = (4.2 ± 0.4) • 10- GeV 
VSvrW^; °M 2 A -M 2 {m 2 * 



where a ~ 1/137 is the electromagnetic coupling constant, and the error has been 
estimated on account of the large- N c approximation of the model. Equations (|73j) fix 

13 See Ref. [22] for a more elaborated version of this method applied to QCD, in which also a p 1 is 
included. 
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the parameters in Eq. (J72J) to be 

M v = 0.70 ± 0.01 GeV , M A = 1.00 ± 0.03 GeV 

F v = 122 ± 6 MeV , F A = 84 ± 7 MeV . (74) 

Feeding expression (J72j) with these values one has a prediction for the full vacuum 
polarization function. In the Euclidean region the prediction is very good for a broad 
range of intermediate values of q 2 . For example, at q 2 = — 1 GeV 2 the error between 
I% 2 ) and U MHA (q 2 ) is only ~ 3% 14 . 

Expanding Eq. (|T2*|) in inverse powers of q 2 , we obtain the following prediction for 
the OPE coefficients: 

a MHA = _^ g ± Q 3) x 1Q -3 Ge y6 ) a MHA = ^ ± Q ?) x 1Q -3 Ge y8 ? ^ 

to be compared to the results in Eq. (|46|) . Again, no dramatic improvement is found as 
compared to previous methods. As with other methods, the prediction for a 8 is worse 
than that for ciq. 

In summary, the errors in the determinations of OPE coefficients from the different 
methods tend to be large in the present model. No method can be claimed to be more 
precise than the others, and the error made is comparable in size to the spread of values 
among the different methods. 

It is difficult to say what happens in the case of QCD. For a 6 , the spread of values 
in QCD found in the literature is not as bad as in the case of our model. This could 
be related to the fact that the size of the coefficient b$ relative to a 6 is much smaller 
in QCD than it is in the model. But also in QCD the discrepancies are clearly more 
serious for a 8 and higher coefficients [T7j . 

One of the lessons from our analysis of the model is that one should not ignore 
duality violations represented by the T> terms in Eqs. f!59H63|) . In Fig. 9, we have 
plotted these T> terms (solid lines). As one can see, they show zeros at certain values 
of s ~ 1.5, 2.4, 3.4, ... GeV. These zeros correspond to values of s where the duality 
violations vanish and they are therefore the optimal points which one would like to 
use as duality points in any finite-energy sum rule analysis. However, the location 
of these zeros is not precisely determined because, as Fig. 9 shows, the zeros are not 
universal but move from one spectral moment to the next. This is one of the reasons 
why the finite-energy sum rule determination of OPE coefficients in Eqs. (j64H68|) has 
large errors. Consequently, as such duality violations are likely to be present in QCD 
as well, we will now use the model to show how these T> terms may be included in a 
realistic analysis also in QCD. 

Our basic assumption will be that our result (p)H|) together with the functional 
form of a damped oscillation as in Eq. (jBTIj) are rather generic results which go beyond 
our particular model. This assumption seems reasonable 15 as these general features 
basically depend on large N c and on the hierarchy of exponentials (J56|) controlling the 
different contributions in Eqs. (|52H55|) . In turn, this depends only on the fact that 
large- N c duality violations are concentrated on the real axis in the complex q 2 plane. 

14 See also Ref. [33] • 
15 See also Ref. 0. 
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Therefore, let us start by assuming a duality violation of the form given in Eq. (J70|) . 
with parameters k, 7, a and (3 to be determined. We remind the reader that in 
Eqs. (I59I60|) the b terms are suppressed by inverse powers of s . Therefore, in a first 
approximation, it may be reasonable to completely neglect the b terms in these two 
equations and, through Eq. (jHHl), extract these parameters k, 7, a and (3 by means of 
a simultaneous fit to Eqs. ()59I60|) . Then the rest of coefficients a6,8,... an d &6,8,... can be 
obtained from Eqs. (j61H63)) in a straightforward way. 

Just to illustrate how this would work in our model, we have determined the param- 
eters in Eq. (fTUj) by doing this simultaneous fit in the window 1.5 GeV 2 < s < 3.5 GeV 2 
with 20 equally-spaced points, with the following result: 

k = 0.026 , 7 = 0.591 GeV" 2 , a = 3.323 and (3 = 3.112 GeV" 2 . (76) 

Fig. El shows in the upper row the result of the fit to Eqs. (|59I60|) . As one can see, the 
overall fit is quite good. And since it improves in the intermediate region of the fitting 
window where both T>^ and X>M have a zero, this suggests that a good strategy is to 
determine the OPE at these zeros. This is what we will do next. 

In the lower row of Fig. El the dashed curve is the result for pt 2 ' 3 l(so) from the fitted 
function (fTUj) and the solid curve is the result for the left-hand side of Eqs. ()61I62|) . as 
calculated from the true values of the model. 

Armed with Eq. (J70|) and the values of the parameters from the fit (J76j) . we predict 
a duality point Sq at which U^(sq) = 0, within the window of our fit, to be at Sq = 
2.350 GeV 2 . Although similar, this value is not quite the same as the one obtained in 
the finite-energy sum rules analysis before, which was 2.363 GeV 2 . This small difference 
has a large impact on the determination of the OPE due to the steepness of the slope in 
the plot of Fig. 9 (lower curves), resulting in the large error quoted in the finite-energy 
sum rule analysis of Eqs. (|67I68|) . 

At this s* = 2.350 GeV 2 , on the other hand, one finds from Eq. fEEJ) that the OPE 
contribution is given by 

a 6 + &alog| + ^-^ = f'dtepit) (77) 

= -0.00251 GeV 6 , 

whereas the exact number for the combination on the left-hand side of Eq. ()77|) obtained 
in the model is —0.00281 GeV 2 , i.e. a ~ 10% error. 16 The same analysis can be repeated 

16 Wc have repeated the analysis including the true values for the coefficients be and bs in the fit to 
Eqs. H59I60(I . In this case the error becomes 14%. 
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Figure 9: Upper row: Plot of the left-hand side (solid line) and fit of the right-hand side 
(dashed line) of Eqs. 15916(J\) . neglecting the b coefficients, as a function of so (in GeV 2 ). 
Lower row: Plot of the left-hand side (solid line) and right-hand side predicted by the fit 
(dashed line) of Eqs. as a function of sq (in GeV 2 j. 



for the combinations appearing in the rest of Eqs. (J621 For instance, we obtain 
a 8 + b 8 \og^-b 6 s* + b -^ (7* 



- dt t 3 p{t) = 0.00329 GeV 8 (s* = 2.307 GeV 2 
Jo 



and 

n + h w s *° ^g 5 + bsS *° 4 bl ° s *° 3 + buS *° 2 h ** + &18 r7Q^ 

a 16 + 16 log— — + — 1 b M s + — (/yj 



o 



- / dt t 7 p(t) = 0.0179 GeV 16 (s* = 2.130 GeV 2 ) 
Jo 



to be compared to the exact values on the left-hand sides of Eqs. (I78|) and (J79|) which 
are 0.00344 GeV 8 and 0.0161 GeV 16 , representing ~ 4% and ~ 11% errors, respectively. 
The gain in precision with respect to previous methods is clear. 



4 Conclusion 

We have developed a model of duality violations that shares many properties with 
QCD, including the fact that resonances in the real world have a non-vanishing width. 
Of course, an analysis of duality violations in a model cannot replace a similar analysis 
in QCD. For example, the precision of the determination of combinations of OPE coef- 
ficients in Eqs. (|77H79|) is not necessarily an indication of what can be achieved in QCD. 
However, we do believe that the lessons from our model analysis are also relevant for 
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QCD. First, the properties of the model that lead to duality violations are also present 
in QCD [7,, notably the fact that the OPE in both is (at best) asymptotic. Second, 
because of the semi-realistic properties of the model, our analysis gives a good indica- 
tion of the nature of systematic errors which result from making the various ansdtze 
which have been proposed in the literature in order to deal with duality violations. 
We think that it should prove useful to attempt an analysis of the duality violations 
represented by the (cf. Eqs. ([59Hb3|) ). using an ansatz similar to that of Eq. (J7UJ), 
also in the QCD case. In addition, we would like to call the reader's attention to the 
additional contributions to the "OPE side" of duality relations, represented by the b 2n 
in Eqs. f!59H63|) . While some of the b 2n m ay be small in the case of QCD, there are 
contributions proportional to positive powers of Sq, which enhance the effects of the 
bin coefficients. To the best of our knowledge, this issue has not been taken into ac- 
count in any of the previous determinations in the literature of higher-dimension OPE 
coefficients from duality. In a realistic analysis of the data, the procedure probably 
should be refined to determine all the a and b coefficients of the OPE, as well as the 
parameters K,j,a and (3 of Eq. (fTOj) parameterizing the duality violations, by means 
of a simultaneous fit to all Eqs. (|59H63|) . in a manner similar to that of Ref. [2D] . 

Finally, in real life one will also have to deal with experimental error bars, of course. 
For example, the experimental error bars in Fig. (jSJ) are large in the large-so region, and 
this will further complicate a precision determination of OPE coefficients. However, 
we hope to have made it clear that much can be learned from considering models such 
as the one analyzed in this paper, and that some of the lessons can be applied to the 
real- world case of QCD. At the very least, an analysis like the one presented here offers 
a way to estimate systematic effects. 
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Appendix 

The OPE of our model stems from the asymptotic expansion of the ip{z) function in 
Eq. ()26|). Using the identity for Bernoulli numbers 

B 2n = 2 C(2n) , n = 1,2,3,... (80) 
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where C{2n) is Riemann's zeta function, one can estimate the minimal error, 5, for this 
expansion as 

2n*z 2n * y ' 

where n* is the optimal value in the sum (|2fij) . This n*, in turn, can be estimated as 
the value at which two consecutive terms in the sum (j2T?j) are approximately equal. 
This results in n* ~ ttz. 

Consequently, after using Stirling's expression for the factorial, one obtains for the 
minimum error in the asymptotic expansion (|2fjjl 

5 ~ e~ 2wz , (82) 

where, here, z ~ ( a^) a ^ l ar S e \q 2 \i an d C = 1 — ^fr f° r finite ^c- The estimate 
(JH2*|) is also valid for the function Ilvv ~ ^aa, because there is no cancelation in the 
difference after inclusion of the right pre-factors. The result (|H2*|) is what appears in 
Eqs. (J2BJ) and (JI2J) in the main text. 

We use the opportunity to also give some details of the derivation of Eq. 1)5 5 j) . From 
its definition (jHSJ), £coLt.( s o) can be written as 



S3) 



The two integrals can be estimated by using the large-g 2 behavior of A(g 2 ), if s is 
large. Referring back to Eq. ()43jl . we note that 



cot 



A 2 " +7r A^ 



2e"^ sin(2x) ±i(l + 2e~ 2y cos(2x)) + O (e~ Ay ) , (84) 



where the plus (minus) sign should be used on the positive (negative) imaginary axis, 
and where we defined 



X = 71 



71 



A 2 



A 2 



C 2 
771 

cos(Cvr/2) + vr— , (85) 

C 

sin(C7r/2) . 



It is then straightforward to estimate each integral by using the bounds | sin(2x)| < 1, 
| cos (2a?) | < 1, and taking £ — ► 1. The unsuppressed ±i in Eq. (|84"j) cancels between 
the vector and axial channels, and we obtain the result given in Eq. (|55p. 
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